**FULL TITLE** 

ASP Conference Series, Vol. **VOLUME**, © **YEAR OF PUBLICATION** 
**NAMES OF EDITORS** 



Shell Galaxies: Dynamical Friction, Gradual Satellite 
Decay and Merger Dating 

Ivana Ebrova, 1,2 Bruno Jungwiert, 2 Gabriela Canalizo, 3 Nicola 
Bennert, 4 and Lucie Jflkova 5 

1 Faculty of Mathematics and Physics, Charles University in Prague, 
Ke Karlovu 3, CZ-121 16 Prague, Czech Republic 

2 Astronomical Institute, Academy of Sciences of the Czech Republic, 
Bocni II 1401/la, CZ-141 31 Prague, Czech Republic 

3 Institute of Geophysics and Planetary Physics & Dept. of Physics, 
University of California, Riverside, CA 92521, USA 

A Dept. of Physics, Univ. of California, Santa Barbara, CA 93106, USA 

5 Dept. of Theoretical Physic and Astronomy, Faculty of Science, 
Masaryk University, Kotldfskd 2, CZ-611 37 Brno, Czech Republic 

Abstract. With the goal to refine modelling of shell galaxies and the use of 
shells to probe the merger history, we develop a new method for implementing 
dynamical friction in test-particle simulations of radial minor mergers. The 
friction is combined with a gradual decay of the dwarf galaxy. The coupling of 
both effects can considerably redistribute positions and luminosities of shells; 
neglecting them can lead to significant errors in attempts to date the merger. 



1 Shells as Probes of the Host Galaxy Merger History 

Shell galaxies contain faint arc-like stellar features. It is widely believed that 
shells are a signature of a merger experienced by the host galaxy. They contain 
at most a few percent of the overall galaxy luminosity, and their contrast is 
usually very low. The model of a radial me rger of a giant elliptical with a 
smal l er galaxy (a spiral or a d warf elliptical) ( Quinn|[l984l ; Dupraz &; Combes] 



1986; iHernquist Quirmlll988l ) seems to be the most successful in reproducing 



regular shell systems. When a small galaxy enters the sphere of influence of 
a giant elliptical on a close-to-radial trajectory, it disintegrates and its stars 
begin to oscillate in the potential of the giant. At their turning points, where 
the stars tend to spend most of their time, they pile up and produce arc-like 
enhancements in the luminosity profile of the host galaxy. 

Attempts to date a merg er from observe d posit ions of shells have been made 
in previous works. Recently, ICanalizo et al.1 ( 20071 ) presented HST/ACS obser- 



vations of spectacular shells in a quasar host galaxy (Fig.Q]) and, by simulating 
the position of the outermost shell by means of restricted N-body simulations, 
attempted to put constraints on the age of the merger. They concluded that 
it occurred a few hundred Myr to ~ 2 Gyr ago, supporting a potential causal 
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Figure 1. Deep HST/ACS images of the host galaxy of t he quasar MC2 
1635+119, so far the only known shell galaxy with a quasar (|Canalizo et al.l 
l2007t iBennert et al.ll2008T ). The left panel shows the original image, the right 
one the residual after the subtraction of the fitted smooth light profile. 



connection between the merger, the post-starburst ages in nuclear stellar popu- 
lations, and the quasar. A typical delay of 1-2.5 Gyr between a merger and the 



onset of quasar activity is suggested by bo th N-body simulations (jSpringel et al 



20051 ) and observations (|Rvan et al.ll2008l ). It might therefore appear reassuring 



to find a similar time lag between the merger event and the quasar ignition in a 
study of an individual spectacular object. However, caution must be exercised 
in estimating merger ages from the location of shells (see below). 



2 Dynamical Friction and Gradual Decay of the Satellite 

While the shell formation, once the dwarf galaxy is disrupted, is basically a test- 
particle phenomenon, the gradual decay of the satellite as well as its braking by 
dynamical friction against the primary can considerably affect the energy dis- 
tribution of oscillating stars, and thus the positions a nd the brightness of shells . 
The dynamical frictio n effect was first poi n ted o ut bv lDupraz &; Combesl ( 19871 ) 



and also discussed bv lHernquist Quinn ( 19881). wh i le the gradual decay, with 



friction neglected, was modelled by Heisler 8z White! ( 19901 ). However, coupling 



of these phenomena was never modelled in much detail. Our goal is to improve 
restricted N-body simulations of shells created in minor mergers by a) including 
dynamical friction, b) improving its implementation by avoiding the use of the 
Chandrasekhar formula, c) coupling it to the gradual decay, d) taking into ac- 
count the present state of knowledge of stellar and dark matter distributions in 
both giant and dwarf ellipticals. A detailed description is beyond the scope of 
this paper. Here, we confine ourselves to a simple example of a radial minor 
merger (Fig. [2]), instructive in showing how an observed shell structure could be 
misinterpreted in terms of the merger time scale (and of the relative pre-merger 
motion) if dynamical friction and gradual decay were neglected. 

In test-particle simulations, the Chandrasekhar formula is commonly used 
to include dynamical friction. Its relative simplicity is made possible, among oth- 
ers, by the oversimplifying assumption of homogeneity of the stellar and dark 
matter distributions. To avoid it, we used the axial symmetry of our merger 
configuration to simplify the integrals over impact parameters and velocity dis- 
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tributions so that they can be solved numerically. The mass of the satellite, 
a key quantity for the efficiency of dynamical friction, is gradually lowered in 
proportion to the mass located beyond its evolving tidal radius. 
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Figure 2. Three snapshots of simulations (3.5, 5 and 7 Gyr after the first 
passage of the satellite, coming from the right, through the center of the 
primary) without (upper row) and with (bottom row) dynamical friction and 
gradual disruption (in the first case, the dwarf instantly disrupts during the 
first passage). Only stars of the dwarf are shown. Each box, centered on the 
primary, shows 300x300 kpc. 



The introduction of dynamical friction and gradual decay dramatically 
changes the appearance of shells as can be seen in histograms of particles' galac- 
tocentric distances (Fig.[3l corresponding to central snapshots of Fig. [2]). While 
the position of the outermost shell is not much affected, its brightness is dras- 
tically lowered. The other shells are shifted and new generations of shells are 
added during each successive passage of the dwarf. Easily inferring the age of 
the c ollision is rendered impossible (as already pointed out bv lDupraz &: Combesl 
1987). The shell systems in Fig.[3l both having the outermost shell at +150 kpc, 
are seen 5 Gyr after the first passage of the two galaxies through each other. 
If we observationally identified the leftmost shell (at —80 kpc in Fig.[3l lower 
panel) as being the outermost one, we would mistakenly estimate the merger 
age to be only ~ 2.5 Gyr. We would also wrongly determine the direction from 
which the dwarf came: assuming the classical picture (based on simulations 
without friction and with instantaneous disruption), the outermost shell would 
be located on the side from which the satellite came, so we would conclude it 
went from the left while the opposite is true. 
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3 Conclusions 

Using even the outermost observed shell to date a merger, and basing on it 
a support for a causal connection between the merger and the quasar, is very 
uncertain. Supposedly, the first formed shell (observed as the outermost one if 
still undissolved and bright enough) is the least affected by dynamical friction 
(since it is formed out of stars released during the first satellite's passage) and 
thus the most reliable for merger dating. In our example, this first shell is very 
weak due to the gradual decay of the satellite. If missed in observations, the 
merger age would be underestimated by ~ 2.5 Gyr; in reality, it is twice as old. 





U 


I), 




) 






1 


I 


J 





-150 -10C -50 50 100 150 



Figure 3. Histograms of galactocentric distances of stars (in kpc) originally 
belonging to the dwarf, at 5 Gyr. Top: Instantaneous disruption, no friction; 
Bottom: Gradual disruption plus friction. Distances are measured from the 
center of the primary, and plotted separately for positions on the side from 
which the satellite came and those on the opposite one (plus/minus sign). 
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